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ABSTRACT 

The objective of this research was to analytically and 
experimentally study the capabilities of adaptive material 
plate actuators for suppressing flutter. Piezoelectrics are 
materials which are characterized by their ability to 
produce voltage when subjected to a mechanical strain. 
The converse piezoelectric effect can be utilized to actuate 
a structure by applying a voltage. For this investigation, 
a two degree of freedom wind tunnel model was designed, 
analyzed and tested. The model consisted of a rigid wing 
and a flexible mount system which permitted translational 
and rotational degrees of freedom. Actuators, made of 
piezoelectric material were afixed to leaf springs of the 
mount system. Command signals, applied to the 
piezoelectric actuators, exerted control over the closed loop 
damping and stiffness properties of the leaf springs, A 
mathematical aeroservoelastic model was constructed 
using finite element methods, laminated plate theory, and 
aeroelastic analysis tools. Plant characteristics were 
determined from this model and verified by open loop 
experimental tests. A flutter suppression control law was 
designed, implemented on a digital control computer and 
tested to conditions 20% above the passive flutter speed of 
the model. The experimental results represent the first 
time that adaptive materials have been used to actively 
suppress flutter. It demonstrates that small, carefully- 
placed actuating plates can be used effectively to control 
aeroelastic response. 


INTRODUCTION 

Flutter, nteraction between the structural dynamic and 
the a f amic characteristics of an aircraft causing 
diverge destructive oscillations of motion, has been 
observed and documented since the era of controlled flight 
began. Historically, passive solutions such as increasing 
structural stiffness, mass balancing or modifying 
geometry have been used to prevent this hazardous 
phenomenon. These approaches result in increased weight 
and cost and decreased performance. During the past 
twenty years, there has been considerable research to 
develop active flutter suppression concepts which use 


conventional aerodynamic control surfaces.^* 2 ’ ^ 

An active control approach eliminates most of the weight 
and performance penalties associated with the passive 
approach and additionally provides flexibility in that the 
control law could be varied with configuration or flight 
condition. Active flutter suppression is not yet a 
common practice in today's commercial or military aircraft 
due to several concerns. Flutter is generally of a 
catastrophic nature, therefore a failure of the system could 
affect flight safety. As a result, system redundancy, 
reliability and maintainability are critical issues. To a 
lesser extent, the control surface authority available to 
maneuver the aircraft with the simultaneous 
implementation of active flutter suppression is also a 
concern. To alleviate these concerns, alternatives to 
utilizing the aerodynamic control surfaces for active flutter 
suppression are being studied. 

The use of secondary controllers made of adaptive material 
is one such concept. There are several classifications of 
adaptive materials; the present study focused exclusively 
on the use of piezoelectric materials. Piezoelectricity is 
the ability of a material to develop an electrical charge 
when subjected to a mechanical strain and vice versa. 
Results available from aeroelastic applications of 
piezoelectric materials are very limited. Static 
aeroelasticity has been the subject of investigations by 
T6 7 81 

Ehlers and Weisshaar. ’ " J They conducted analytical 
studies on laminated composite wings with embedded 
actuators, looking at pure torsional and bending 
deformations. They reported that through feedback to 
embedded adaptive material layers, the divergence speed is 
altered, implying also that lift effectiveness is influenced. 
The augmentation or replacement of conventional 
aerodynamic control surfaces with strain actuation for 
aeroelastic control has been the focus of an analytical 

investigation of a typical section by Lazarus, Crawley and 
[91 

Lin. They found that strain actuation via piezoelectric 
elements may provide a viable and effective alternative to 
articulated control surfaces for controlling aeroelastic 
response. Investigation of flutter suppression for lifting 

surfaces and panels has been done by Scott. This 
analytical study considered controlling flutter at supersonic 
speeds using full state feedback. 

The purpose of the present study was to investigate flutter 
suppression using piezoelectric plates as actuators. 
Results from analyses and experiments demonstrating this 
technology are presented. This paper includes: discussion 
of piezoelectric materials; details concerning the 
experimental model, the wind tunnel facilities and the 
testing procedures; details concerning the analytical model 
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development, the analysis results and the control law 
design; and the experimental results and comparison with 
analytical predictions. A more in-depth presentation of 
the material contained in this study may be found in 
reference 11 . 

EIE ZQELECTRIC MATERIALS 

A material which, when subjected to a mechanical load, 
accumulates an electrical charge is said to have 
piezoelectric properties. Many naturally-occuring crystals 
have piezoelectric capabilities but their non-uniform 
properties encouraged research into manufacturing 
materials which would produce electromechanical 
coupling. Certain polymers and ceramics consist of 
crystalline subdomains which are bipolar in nature. 
Piezoelectric properties can be induced in these materials 
by applying a large electrical field across them. This 
induces an orientation of the dipoles such that the positive 
and negative poles of the individual domains are aligned 
with the applied field, denoted the 3-direction, (figure 1 ). 
The orientation remains after the inducing field is 
removed. In order to subsequently deform the material, a 
smaller voltage is applied through electrodes on opposite 
faces of the material. The most common configuration is 
to place the electrodes on the faces parallel to the poled 
axis and to apply the voltage in the same direction as the 
original inducing field. As shown in figure 2, the 
material deforms both through the thickness, denoted the 
^33 effect, ( voltage in the 3-direction and displacement in 

the 3 -direction), and in the in-plane directions, denoted the 
<*31» ( voltage in the 3-direction and displacement in the 1 - 
direction). Applying a voltage through the thickness in 
one direction induces in-plane expansion; applying it in 
the other induces in-plane contraction. 

Piezoelectric plates can be configured in different ways to 
accentuate the displacements or forces being generated. 
The in-plane expansion and contraction of adaptive 
materials may be utilized by bonding actuating plates to 
either side of a center shim, (figure 3). One is expanded 
while the other is contracted resulting in a bending 
displacement much greater than the length deformation of 
either of the two layers. This configuration, which takes 

advantage of the Poisson-like d 31 effect, is referred to as a 
bimorph or a bender element. It will serve as the primary 
actuator mechanization for the investigation described 
herein. 

Lead Zirconate Titanate (PZT), a piezoceramic, was 
chosen for this investigation. The material properties for 
PZT are given in table 1. The material properties have 
been treated in this work with linear relationships. This 
assumption is valid for low applied voltages and small 
deformations, Nonlinearities and nonidealities of these 
materials have been well-documented by references 6 , 7 
and 9. 


EXPERIMENT A APPAPAT^ 

WIND TUNNEE 

The piezoelectric flutter suppression model was tested in 
the Langley Research Center Flutter Research and 
Experiment Device (FRED). The FRED is an open 
circuit table top wind tunnel with a maximum operating 
velocity of 85 miles per hour ( approximately 1500 inches 
per second ). The test section is six inches by six inches, 
and is constructed of plexiglass for model viewing. The 
air is pulled through the tunnel by a 2 horsepower motor 
and smoothed by a single honeycomb screen at the 
beginning of the contraction duct. Models are mounted 
from the removable ceiling of the test section. Figure 4 
shows the test article suspended from the test section 
ceiling. 


The piezoelectric flutter suppression wind tunnel model 
consists of four integral components: a flexible mount 
system, a rigid wing, piezoelectric plate actuators and a 
strain gage bridge. There were three driving factors in the 
test article design: the model had to flutter within the 
wind tunnel envelope; had to fit inside the wind tunnel 
with certain margins of safety; and had to have flat 
surfaces on which piezoelectric actuating plates could be 
mounted. Therefore, a rigid wing was designed with a 
flexible mount system that would reside outside of the 
wind tunnel. The flexible mount system provided two 
degrees of freedom: plunge ( out-of-plane translation) and 
pitch ( rotation about the wing pitch axis). With the rigid 
wing attached to the flexible mount system, the two 
degrees of freedom became coupled by virtue of the wing 
mass distribution. 

The mount system suspended the wing by two pins 
through slots in the test section ceiling and provided the 
plunging and pitching freedoms by virtue of separate 
spring tine mechanisms. The plunge mechanism consists 
of two spring steel plates or tines separated by .75 inches 
and clamped at both ends. The pitch mechanism is a 
single spring tine connected to the wing at the leading 
edge and at the .2353 chord location, wnere there is a 
bearing-like mechanism which allows for free rotation. 

The wing consists of three sections. The primary wing 
structure is formed from one eighth inch thick aluminum. 
It has a chord of 2 inches, with the pitch pivot at the 
midchord. A balsa wood extension overlays the aft half of 
the primary structure and extends the chord length to 4.25 
inches. The trailing edge of this section was coated with 
aluminum to provide a mass ballast. All three sections 
extend the full span of the wing, which is 4 inches. 
Model characteristics are presented in table 2. 

Piezoelectric ceramic plates were installed near the root of 
one of the steel plunge spring tines to actuate the test 
article. Two plates, 1-1/2 inches long and 1 inch wide, 
were bonded to opposite sides of the plunge spring tine. 
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with their poles both oriented towards the steel, to form 
an actuator. The .0075 inch thick plates were electrically 
isolated from the steel by the bonding layers. Small 
copper tabs afixed beneath the plates during the bonding 
process served as the means of applying voltages to the 
bonded-side electrodes, (figure 5). 

A strain gage bridge was mounted near the base of one of 
the spring tines, with two gages on either side of the tine. 
The gages were configured to measure the cantilever 
bending strain and were used as the feedback signal for the 
controller. 

DIGITAL CONTROTTirq 

The control laws were implemented using a personal 
computer, with a 80386 processor and 80387 co-processor 
running a real time Unix operating system. The control 
laws are programmed in the C-language and use floating 
point arithmetic for all control law calculations. The data 
acquisition system uses 12 bit analog-to-digital 
converters. Additional details of the controller and data 
acquisition system can be found in Reference 12. 

ADDITIONAL INSTRUMENTS^ 

Aside from the strain gage, there were two sensor systems 
m use: a linear accelerometer and a hot wire anemometer. 
The accelerometer was used in system identification: 
experiments; it served as a roving measurement, being 
placed where applicable for different experiments. During 
zero airspeed testing, it was located on the wing, however, 
during flutter testing it was installed on the clamping 
block of the mount system. A Kurz 443M air velocity 
meter, a hot wire anemometer, provided visual readouts of 
the test section airspeed. The probe was inserted into the 
flow just behind the model in the test section. Airspeed 
was measured by removing the model from the tunnel to 
eliminate blockage and the influence of wing oscillations 
on the reading. 

An Apex Microtechnology P83A operational amplifier 
was used to boost the input voltage to the piezoelectric 
actuators. The signal source or input voltage was 
amplified by a factor of 25, with a limit on the output 
voltage equal to the power supply voltage, which in this 
experiment was +/- 80 volts. 

ANALYSES 

The analytical model was developed incorporating the 
aeroservoelastic equations of motion, a model of the 
control computer dynamics and experimentally-determined 
correction factors. Analyses were performed utilizing the 
aeroservoelastic equations of motion. Results from 
system identification tests, which include natural 
frequencies, structural dampings, and scale factors for the 
actuator and the sensors, were incorporated into the 
aeroservoelastic equations of motion prior to the open 
loop analysis. Aeroservoelastic equations of motion based 
on Lagrange’s energy equations represent a summation of 


forces which include the inertial, dissipation, and internal 
restoring forces, and the reduced-frequency dependent 
aerodynamic forces due to the structural motions. 
Classical techniques available to the aeroelastic 
community were utilized to develop the equations. 

MODEL DRSIfiN 

Structural Mnrf«»| A finite element model of the wind 
tunnel wing and the two degree of freedom mount system, 
(figure 6), was constructed and analyzed. The wing was 
modeled with solid elements and concentrated masses. 
The spring tines were modeled with plate elements. A 
torsional spring was added at the pivot point to better 
represent the experimentally-determined pitch frequency. 
The primary purpose of the finite element model was to 
generate the structural mass and stiffness matrices. The 
model was also utilized to perform parametric studies for 
the design of the mass ballast and to determine the 
placement of the actuating plates necessary to obtain the 
maximum control effect. 

A normal mode analysis was performed to generate the 
two important natural frequencies and mode shapes of the 
model. The first mode, designated plunge due to the 
dominance of translational motion, was predicted at a 
frequency of 7.8 Hz. The second mode, which is 
characterized by the pitching of the wing relative to the 
mount system, has a natural frequency of 10.9 Hz. 

Aerodynamic Model Unsteady aerodynamics were 

calculated using the Doublet-Lattice Method [13 ^ as 
implemented in the Aeroelastic Vehicle Analysis (AVA) 
system of computer codes. AVA uses the modal 
displacement vectors to calculate the generalized 
aerodynamic forces (GAFs) at discrete reduced frequencies. 
The program output is a table for each reduced frequency, 
where the columns of the table correspond to modal 
deflections, and the rows correspond to modal pressures or 
forces. The aerodynamic model has 5 chordwise boxes and 
10 span wise boxes for a total of 50. The GAFs were 
calculated at Mach .05 for 8 values of reduced frequency 
ranging from .001 to 2.0. 

QPEN LOO P AEROELA STIC analysis 

The structural model and the aerodynamic model were 
combined to form the open loop aeroelastic equations of 
motion. Flutter analyses were conducted by analyzing 
this model for sea level density using the STABCAR 

[Hi 

computer code. The velocity root locus plot of figure 
7 shows the open loop flutter characteristics. The plot 
traces the eigenvalues of the system as the airspeed is 
increased. The imaginary axis represents the point of 
neutral stability or zero damping, where, theoretically, 
responses will neither converge nor diverge. Flutter, 
defined as an oscillatory divergence, is represented on a 
root locus plot by an eigenvalue crossing this axis into 
the right half plane. 
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Figure 7 shows the coaJescent behavior of the plunge and 
pitch inodes for sea level density as the airspeed increases. 
The frequencies of the two modes migrate towards one 
aether as the aerodynamics couple the their motion. 
When the frequencies are close together, the modes interact 
with one another and the system is driven unstable, shown 
7 ^ tower frequency mode eigenvalue crossing into the 
right half plane. Flutter is predicted to occur at a velocity 
of 560 inches per second and at a frequency of 9. 1 Hz. 


having only displacement measurements/ 18 ^ The cunent 
control law utilizes these dynamics of the implementation 
scheme, requiring only gain feedback. Figure 9 shows 
that for a 20 Hz sample rate, the phase is -270 degrees, or 
+90 degrees at the Nyquist frequency, 10 Hz. Thus the 
phase characteristics simulate a derivative near the 
anticipated flutter frequency, which is the frequency of 
interest for control. 1 


STATE SPACE A E ROSFRVQFT.ASTTr C ONTROL- LAW DESIGN CLOSED LOOP 

EQUATIONS QF MOTION A NALYSIS 


The formulation of a state space aeroservoelastic model 
requires that the piezoelectric control forces be combined 
with the aeroelastic equations of motion described above. 

The control forces from the piezoelectric elements were 
developed utilizing laminated plate theory in conjunction 
with Lagrange s energy method. The modeling of 
piezoelec trical systems requires consideration of both 
mechanical and electrical behavior. Coupling between 
mechanical stresses and electrical fields is analytically 
represented by constitutive relationships which contain 
both the electrical quantities and the mechanical quantities. 
This relationship is utilized in expressing the potential 
energy- the presence of the strain due to electrical energy 
generates a term which is not traditionally observed in the 
mechanical potential energy expression. The control force 
is calculated using a finite difference program. 

The aerodynamic forces within the aeroelastic equations of 
motion produced by the Doublet-Lattice code are 
transcendental functions of reduced frequency. To 
incorporate them into the state space equations of motion, 
they must be approximated by rational functions of the 

Laplace variable, s/ 15] The Integration of Structures, 

Aerodynamics and Controls ( ISAC ) system of codes 
was used to generate the resultant s-plane GAFs. 

The development of the first order state equations of 
motion was accomplished using MatrixX, a 
commercially-available software package from Integrated 

Systems Incorporated/ ^ These continuous system 
matrices were then discretized using a 20 Hz sample rate. 

A block diagram of generic form, describing the closed 
loop system is presented in figure 8. 

MODELING THE CONT ROL COMPUTED 
DYNAMICS 

Illustrated in figure 9 for again feedback control law, the " 
control computer introduces its own dynamics into the 
feedback path. The digital controller implementation 
scheme shifts the output data by one sample and applies a 
zero order hold. The frequency response of the digital 
controller is different for different sample rates. By using 
a sample rate which provides the correct amount of phase 
shift near the frequency of control interest allows the 
system to simulate analog derivative feedback, despite 


The control law design requirement for the current study 
was to suppress flutter for the largest velocity possible. 
Due to the simplicity of both the design objective and the 
test article, the aeroelastic phenomenon was controllable 
through a smgle input / single output control law. Gain 

feedback utilizing the dynamics of the discretization 
process was employed. The strain-proportional voltage 
was the input to the control law. The implementable 
feedback gain was limited to +/- 33 due to a voltage 
limitation of the amplifier driving the piezoelectric plates. 

Control law design is traditionally performed in the 
continuous domain, however because the control law 
computer discretization was an integral part of this design, 
the discrete domain model was utilized in this 
investigation. On a diagram of discrete system 
eigenvalues, the stability condition corresponds to the 
location of the roots relative to the unit circle. Roots 
located outside of the unit circle correspond to 
instabilities. That is, the imaginary axis of the 
continuous complex plane, maps to the unit circle of the 
discrete complex plane. 

Design models woe constructed from the aeroservoelastic 
equations of motion at several velocities. The model 
representing the system at 580 inches per second, a 
condition above the flutter prediction, was the initial 
design model. A discrete gain root locus, constructed by 
varying the gain from 0 to 120 is shown in figure 10. 
Each eigenvalue trace of the discretized system begins at 
the open loop system values which correspond to a 
feedback gain of zero. One pair of roots shown in figure 
10 is unstable for the open loop case as this velocity 
corresponds to a condition above open loop flutter. The 
flutter mode eigenvalues stabilize for small feedback 
gains; they migrate inside the unit circle almost 
immediately. As the gain increased, the eigenvalues 
continue to migrate in the unit circle for feedback gains up 
to 108, where one destabilized again. 

the stability criterion can be expressed as a limit on the r 
magnitude of the eigenvalues. The magnitude of the 
largest eigenvalue must be less than 1.0 for the system to 
be stable. Figure 11 shows the value of the maximum 
magnitude of the eigenvalues plotted against feedback 
gain. The design model, linearized at 580 inches jte? 
second, is shown to be stable for gains between 14 and 
108. This model, however, represents the system at only 
one airspeed. The same figure shows the variation with 
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gain for several airspeeds. For 1 inch per second airflow, 
the model was stable for all gains below 104. The 
maximum airspeed for which gain feedback stabilized the 
system was found to be 1300 inches per second. As 
indicated in the figure, gains of between 103 and 108 
stabilized the system at this airspeed. If physically 
attainable and no stability margins are required, 103 is the 
optimal gain. Practical limitations on the gain as 
described previously, however, did not allow a gain of this 
magnitude to be implemented. 

Using the largest allowable feedback gain, 33, a velocity 
root locus was constructed. Figure 12 shows the traces of 
the eigenvalues as the velocity is increased from 0 to 700. 
As with the gain root locus, a stable system has 
eigenvalues all lying within the unit circle. The figure 
shows the potential flutter mode migrating towards the 
unit circle for increasing velocities. The root crosses the 
stability boundary, predicting closed loop flutter at 648 
inches per second. The previous graph, (figure 1 1), shows 
the variation of maximum eigenvalue magnitude as a 
function of gain for this velocity. The influence of 
increasing velocity can also be seen on this graph by 
examining the dashed vertical line representing a feedback 
gain of 33. For a velocity of 648 inches per second, the 
eigenvalue trace intersects the stability boundary. Recall 
that this is the predicted closed loop flutter speed. For 
velocities greater than the closed loop flutter speed, there 
is a substantially increased magnitude of instability. This 
is illustrated by the data corresponding to a velocity of 
1300 inches per second. 

EXPERIMENTS 

Experimental results from bench tests and wind tunnel 
tests are presented in this section. Bench tests provided 
zero airspeed system identification results. The wind 
tunnel results to be discussed include system identification 
and flutter testing. Open and closed loop flutter tests were 
conducted and the results are compared to one another as 
well as to analytical predictions. 

SYSTEM IDENTIFICATI ON TESTING 

System identification testing was performed using several 
techniques to extract modal frequencies, dampings, transfer 
functions and general system behavior. Impulse response 
functions of the accelerometer were generated by hammer 
taps. A second technique, employed to obtain a more 
dramatic response from the pitch mode was to pluck the 
spring tine and record the free decay data. 

'• -T” t : 7:VV ??'•'='- • - 

Impulse Tests Hammer taps to produce impulsive 
acceleration responses were used at zero airspeed to extract 
the undamped natural frequencies. Figure 13 shows time 
histories of hammer input and accelerometer response and 
the normalized power spectral density of the acceleration 
response due to the hammer input. The power spectral 
density indicates that the natural frequencies of the final 
configuration are 7.9 and 11.1 Hz. The structural 
damping of the plunge mode was also determined from 


this data by taking the ratio of the frequency width of the 
peak at the half amplitude and the natural frequency. The 
damping ratio is half of this value, .017. 

This identification method could be used only at zero 
airspeed because the wing could not be directly accessed 
while in the wind tunnel. Additionally, the amount of 
disturbance introduced to the flow by the presence of the 
accelerometer and its lead wire drastically altered the 
aerodynamic behavior. 

Free Decay Tests Free decay tests were used to extract 
the open loop damping of the pitch mode and to compare 
the open and closed loop dampings of the plunge mode. 
This free decay technique was effective at low airspeeds, 
where the plucking did not perturb the model enough to 
induce large oscillations. 

The first set of tests were conducted by constraining the 
model in plunge and then plucking the pitch spring tine. 
This method was utilized because the impulse response 
testing failed to extract data suitable for determining the 
pitch mode damping. Random input tests could not be 
used as there was no actuator for this degree of freedom. 
No transfer functions between the input and output can be 
derived using this method because the input signal can not 
be recorded. The usefulness of this test was in the 
determination of the pitch mode damping and frequency. 
A resultant time history from the free decay testing is 
presented in figure 14. Acceleration response of the open 
loop system was generated by applying 7 impulses during 
18 seconds. Each response was decayed before the next 
was applied. This data was analyzed using the logarithmic 

decrement technique. The pitch mode damping ratio, C* 
which is half of the structural damping, g, was determined 
to be .055. 

A similar test was performed to compare the open and 
closed loop behavior at zero airspeed without the plunge 
mode being constrained. The strain response to pluck 
tests of the open and closed loop systems were compared 
in figure 15. Both responses were normalized such that 
the magnitude of the first peak was 1.0; the data were 
obtained using 20 Hz sample rates. From this plot, the 
damping is shown to have been increased by the presence 
of the controller. 

OPEN LOOP FLUTTER TFSTTNO 

Flutter testing was conducted by placing the model in the 
wind tunnel at zero airspeed and then increasing the 
velocity. Data was taken for several minutes at various 
airspeeds. During this time interval the model sat at the 
tunnel condition. Flutter points were defined as the 
lowest airspeed at which the magnitude of the oscillations 
diverged within this time interval. The turbulence within 
the tunnel was relied upon to be sufficient to perturb the 
model. The increments in velocity were made smaller and 
data taken more frequently as the speed got near the 
predicted flutter value. Flutter was encountered 
experimentally at 580 inches per second. A time history 
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of the strain gage during a run in which flutter was 
encountered, (figure 16), shows the divergent oscillations 
which begin growing at 4.5 seconds and continue to grow 
until the maximum possible amplitude was reached at 9.0 
seconds. At this amplitude, safety stops of the tunnel 
inhibited the model's motion so that it wouldn’t be 
destroyed. Frequency domain analysis, performed by 
taking fast Fourier transforms of the strain time history 
data indicates the flutter frequency of 9.4 Hz. 

CLOSED LOOP TF.STTNC 

Closed loop flutter testing was initiated by activating the 
control law at zero airspeed. Proceeding in the same 
manner as the open loop flutter testing, the speed was 
increased until flutter was encountered. 

The closed loop data were compared with open loop data. 
Figure 17 shows time histories of the strain response for 
the two cases. These data were obtained just below the 
open loop flutter speed, at approximately 575 inches per 
second. The magnitude of the peaks for the closed loop 
case are decreased below the magnitudes of the open loop 
peaks. Due to limitations in the controller programming 
the controller update rate and data sampling must be 
consistent. The closed loop data was therefore obtained 
using a 20 Hz sample rate. This low sample rate is 
responsible for the spiky look of the graph. 

The results of the flutter experiments and analyses are 
summarized in figure 18. The top graph of figure 18 
shows the natural frequency variation as velocity is 
increased. Examining the zero airspeed data, the measured 
open loop values (7.9 and 11.1 Hz) are slightly higher 
than the analytical predictions of these frequencies (7.8 and 
10.9 Hz). The higher frequency mode is the pitch mode; 
the lower frequency mode is the plunge mode. The 
analytical values for the two frequencies approach each 
other as the velocity increases. When frequencies coalesce, 
flutter generally onsets soon thereafter. The experimental 
open loop flutter frequency, 9.4 Hz, is slightly larger than 
the frequency which the analysis predicted for flutter. The 
final point on this graph is the experimental closed loop 
flutter frequency, 9.7 Hz. 

The lower graph of figure 18 shows the damping ratio 
variation with velocity. Flutter is indicated on this graph 
by the damping ratio of one of the modes going to zero. 
Examining the zero airspeed data, the measured values of 
damping ratio (.017 and .055) were utilized in the 
analytical model. The analysis results are plotted over the 
velocity range until flutter occurs, (i.e. the lower mode 
crosses the zero damping axis), at 560 in/sec. This value 
is 3.5% conservative as the experimental flutter speed was 
580 in/sec. The closed loop data at zero airspeed indicates 
that the control law provided an increase in the damping. 
Only the data for the plunge mode is presented. 
Analytical predictions of the damping behavior are 
presented over the entire velocity range, until closed loop 
flutter is predicted to occur at 648 in/sec. The 
experimental value was 697 in/sec. The predicted value is 
7% conservative. Using single input/single output gain 


feedback, the flutter speed was predicted to improve by 
15.7%. The controller actually achieved an improvement 
of 20% of the flutter velocity. 

CONCLUSIONS & RECOMME NDATIONS 

This research effort has resulted in the first experimental 
demonstration of flutter suppression employing 
piezoelectric actuators. A digital control law was designed 
and implemented based on a discretized model. Open and 
closed loop flutter tests were conducted, with excellent 
correlation achieved by analytical predictions. The open 
loop flutter velocity prediction was 3.5% conservative 
while the closed loop prediction was 7.6% conservative. 
The analysis indicated that the flutter velocity would be 
increased by 15.7% with the control law; a 20% 
improvement was observed experimentally 

It is recommended that further research be performed in the 
area of controlling the aeroelastic responses of a vehicle 
utilizing piezoelectric actuators. General areas which have 
yet to be tackled include buffet, gust, and maneuver load 
alleviation systems. Further vibration suppression 
research needs to be conducted using a more realistic and 
complex model which incorporates strain-actuating 
elements within a flexible wing design. The concept has 
been proven to work, however, it has not yet been shown 
to be workable in terms of real aircraft Experiments on a 
larger scale are now called for. 
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Property 

Value 

d 31 

166 x 10'* 2 m / V 

d 33 

-360 x 10' 12 m / V 

Maximum Voltage 

.7 x 10 6 V / m 

Young's Modulus 

63 x 10 9 N / m 2 

Poisson's Ratio 

.33 

Density 

7.65 x 10^ kg / 


Table 1 Materia! Properties for Lead Zirconate 
Titanate (PZT) 


PROPERTY 

VALUE 

Wing Mass 

.09 Ibra 

Wing I a 

.134 lbm-in2 

Wing Pitch Axis 

25. % chord 

Wing Center of Gravity 

36.6 % chord 

Plunge Natural Frequency 

7.9 Hz 

Pitch Natural Frequency 

11.1 Hz 

Plunge Mode Structural Damping 

.034 

Pitch Mode Structural Damping 

.11 


Table 2 Measured Model Characteristics 
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Fig. 1 Electrode placement on piezoelectric plate. 


d 31 «ffoct 



Fig. 2 Thickening (d^) and foreshortening (d 3 2 ) effects 
of applied voltage. 


Fig. 4 Photograph of the test section with the wind-tunnel 
test article suspended from the ceiling. 
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Fig. 3 Bimorph or Bender configuration of piezoelectric 
plate elements. 


•v 


Fig. 5 Piezoelectric actuator attachment to plunge 
spring tine. 
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Fig. 6 Finite element model. 
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Fig. 9 Frequency response plot of implicit controller 
dynamics for 20 Hz sample rate (time delay and 
zero order hold). 


Fig. 7 Root locus-eigenvalue locations as velocity 
is varied. 



Fig. 8 Block diagram of closed-loop system. 


Imaginary Part 



Fig. 10 Gain root locus of the discrete system at open 
loop flutter, 580 inches/sec with gains 
from 0 to 120. 
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Fig. 1 1 Maximum complex absolute value as a function 
of feedback gain for various airspeeds. 
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Fig. 13 Results of impulsive input system identification 
tests. 


Real Part 
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Fig. 12 Discrete system root locus for varying velocity Fig. 14 Time history of acceleration response to plucks 
of the closed loop system with feedback gain at zero airspeed (model constrained in plunge). 








NORMALIZED 
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Fig. 15 Open and closed loop strain response to free 
decay test at zero airspeed 
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Fig. 16 Open loop flutter point, 580 inches per second. 


Fig. 18 Analytical and Experimental Flutter Results. 




Fig. 17 Experimental open and closed loop strain 

response to wind tunnel turbulence just below 
flutter velocity (575 inches per second). 
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